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Abstract 

We obtain a fundamental measure density functional for mixtures of parallel hard cylinders. To this 
purpose we first generalize to multicomponent mixtures the fundamental measure functional proposed by 
Tarazona and Rosenfeld for a one-component hard disk fluid, through a method alternative to the cavity 
formalism of these authors. We show the equivalence of both methods when applied to two-dimensional 
fluids. The density functional so obtained reduces to the exact density functional for one-dimensional 
mixtures of hard rods when applied to one-dimensional profiles. In a second step we apply an idea put 
forward some time ago by two of us, based again on a dimensional reduction of the system, and derive a 
density functional for mixtures of parallel hard cylinders. We explore some features of this functional by 
determining the fluid-fluid demixing spinodals for a binary mixture of cylinders with the same volume, and 
by calculating the direct correlation functions. 

PACS numbers: 61.20.Gy,61.30.Cz,64.75.+g, 
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I. INTRODUCTION 



The fundamenta 
spheres (HS) 



measure density functional originally derived by Rosenfeld for a fluid of hard 
50 can be considered as the most sophisticated density functional (DF) 



that has been successfully applied to the study of the highly confined HS fluid [a, 



70 and to the 



HS freezing [81. The theoretical formalism initially developed by Rosenfeld to obtain the HS fun- 

nnn 

damental measure functional (FMF) HI, |2J, |3fl was later complemented with the concept of dimen- 
sional crossover to zero dimension to obtain a DF that adequately describes the HS freezing ^.[loll. 
By dimensional crossover it is understood that when reducing the dimension of a D-dimensional 
system to that of a D'-dimensional one (e.g. by confining), the D-dimensional DF crosses over to 
the D'-dimensional one. This zero-dimensional (OD) crossover was later employed by Tarazona 
and Rosenfeld to introduce a new cavity formalism which showed how this unique property, to- 
gether with the exact expressions for the zero- and one-dimensional (ID) HS functionals, are the 
only requirements needed to derive their final versions of the FMFs for the hard disk (HD) and the 
HS fluids llHO. These versions which, as it is standard in all FMFs, assume that the density profile 
dependence enters in the functional only through a finite set of weighted densities, leave little free- 
dom for improvements without destroying the important dimensional crossover property. Recent 
efforts have been made to derive a HS functional with an imposed equation of state (EOS) as its 



uniform limit 11121 . 



13|]. This imposed EOS (for instance the Carnahan-Starling EOS ) describes the 



fluid better than the scaled particle result [14] (the uniform limit of all the FMF derived from first 
principles) in the description of the HS liquid. However, all these modifications can be done at the 
expense of losing some, or all, of the dimensional crossovers — a crucial property if one wants to 
study highly confined fluids. 

Thus, we can say that the fundamental measure theory (FMT) is close to its edge in the sense 
that it is questionable that any improvement can be achieved without rendering it intractable II 1511 . 
For instance, it has been argued that the inclusion of an infinite set of weighted densities can 
remove some defects of the HS FMF, because this is what happens — with the addition of a few 



18 



more weighted densities — for the FMF for parallel hard hexagons [|16fl (which is constructed from 
the corresponding functional for parallel hard cubes [|l VL 



19fl). After all, a circle is a polygon 



with infinite 



constructed [20, 



y many sides. Similar conclusions are reached when FMFs for lattice models are 
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220 . In fact, the concept of cavity can be generalized in those models so as to 



account for higher orders of accuracy in the correlation functions 



240- 
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The first extension of the FMF to general anisotropic particles was carried out by Rosen- 



feld Il25ll . However this extension only works for isotropic fluids (with particle axes randomly 
oriented), because the Mayer function is not recovered from the low-density expansion of the 
proposed FMF when particles have a preferred alignment [26]. In recent works the definition 
of the one-particle weights necessary to calculate the weighted densities has been extended in 
such a way as to include an effective dependence on the orientations of both interacting particles 



L2' 



3, 28, 29, 3]. 



This way the Mayer function can be exactly decomposed as a sum of convolu- 



tions between those extended weights. Different FMFs for freely rotating particles: for a mixture 



of HS and hard needles [|27L 12811 . an d for a mixture of hard needles and infinitely thin hard plates 
I30I have been proposed using these extended weights, and the resulting FMFs have some of 
the desired dimensional crossovers. However, the requirement of taking the breadth of the particles 
vanishing small seems to be indispensable to construct such functionals. 

The parallel particle alignment (or the restricted orientation approximation) is a fundamental 
restriction that has to be taken if we want to derive a FMF from first principles, without any 
approximation about the particle characteristic lengths. The FMF for parallel hard cubes [17] and 
its later extension to hard parallelepipeds with orientation restricted to three perpendicular axes 



(Zwanzig model) 1118 



were the first examples of FMFs for anisotropic particles derived form 
first principles using the original Rosenfeld's formalism and adding the dimensional cross-over 
constraint. The FMF for the Zwanzig model was applied to the calculation of phase diagrams 
of the one-component fluids made of hard rods and hard plates fl3 1(1 . and of phase diagrams of 



poly disperse rod-plate mixture Q32J, 133Q. Also, it was applied to the study of interfacial phenomena 
in liquid crystals in three [3^] and two [ 35 ] dimensions. Reference [36] summarizes most of the 
works done on the study of wall-liquid crystal fluid interfaces using the Zwanzig model. 

In this article we want to derive a FMF for another particle geometry with the parallel alignment 
restriction. In this case the fluid is composed of a mixture of parallel hard cylinders (PHCL) with 
different radii Ri and lengths Li. To achieve this we will first of all extend the FMF obtained 
by Tarazona and Rosenfeld [11] for a one-component HD fluid to a multicomponent mixture, 
using an alternative approach to the cavity formalism used by these authors. We prove that the 
resulting functional conforms all the dimensional crossovers and thus we show that both methods 
are completely equivalent. In a second step we derive a functional for a mixture of parallel hard 
cylinders starting from the already obtained HD functional by applying a differential operator, as 
explained in Ref. (|l8l] . This procedure guarantees the dimensional crossover from 3D to 2D (|l8l] . 
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The parallel alignment restriction is, of course, a hard constraint that prevents the use of the 
derived functional in the study of those phenomena governed by changes in the orientational or- 
dering of the constituent particles, as it often occurs in liquid-crystals. However, as the spatial 
correlations are accurately treated, those phases with a high degree of orientation, such as the 
smectic, columnar or crystalline phases at high pressures, should be well described by the present 
functional. Also, the study of non-uniform polydisperse liquid-crystalline phases, which are fre- 



quently present in experiments on colloidal mixtures 11371 13811 . is such a difficult task that the 
parallel alignment simplification seems to be the only way to take some steps forward in that 
direction. 

The article is organized as follows. In Sec. [U we derive the FMF for mixtures of HD. In Sec. Hill 
we use this result to construct a FMF for mixtures of PHCL. Section [IV] discusses the uniform 
mixture of PHCL and in Sec. [V] we derive the direct correlation functions for this mixture. After a 
section summarizing the results of the paper, we include two appendices. In Appendix IA1 we proof 
that the FMF for a mixture of HD has an exact 2D— >1D crossover. This is, of course, inherited 
by the FMF for PHCL. In Appendix [B] the expressions for the geometric terms defining the direct 
correlation functions are explicitly displayed. 

H. FUNDAMENTAL MEASURE DENSITY FUNCTIONAL FOR A MIXTURE OF HARD DISKS 

In this section we will derive a density functional for mixtures of HD based upon Tarazona and 
Rosenfeld's proposal for a one-component HD fluid obtained by using the OD cavity formalism of 



FMT fl 1 IN . We will maintain the functional structure of the excess part of the free energy density, 
$(2D) ^ ^ an( j ex t enc i it t a multicomponent mixture by calculating the kernel Kij (r ) (i and j label 
disk species) which enters the definition of the two-particle weighted density N(r) (see below). 
This functional structure is 

$(2D) (r) = _ no(r) i n [i _ n2 ( r)] + N W , (1) 

l — n 2 (jj 

and its extension to mixtures amounts to writing 

N(r) = J dri J c/r 2 pi(ri)p j (r 2 )%(r - ri, r - r 2 ), (2) 



i,0 



%(n,r 2 ) = ^ 0) (riVf (r 2 )i^-(r 12 ). (3) 
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The one-particle weighted densities n (r) and n 2 (r' 



those of Rosenfeld's for a FMF for HD mixtures lliSyD, i.e 



n 2 (r) being the local packing fraction] are 



«o(r) = E/ ^(r')^ 0) (r-r'), (4) 



n 2 (r) = ^2 J dr'pi(r')u. 



(2) 'r-r'), (5) 



with Pi(r) the density profile of species i and cj^ (r) the one-particle weights defined as 

(o), . 8(Ri-r) ( 2 ), N , . 
w ; ( r ) = —5-^ — > ^ ( r ) = - r ), (6) 

Ri (i — 1, 2, • - • , c, with c the number of components of the mixture) being the particle radii, 
and S(x) and Q(x) the Dirac delta and Heaviside step functions respectively. Equation (0) is the 
natural extension to multicomponent mixtures of the two-particle weighted density N(r) initially 
introduced for a one-component fluid in Ref. Ill ill . The authors of this work found the expression 
for the kernel K(r 12 ) through the requirement that inserting 0D density profiles in the excess part 
of free energy 

/^l 2D) [p(r)] = jdr^\r) (7) 

should recover the interaction part of the free energy of a 0D cavity $ (° D ) = 7V+ ( 1 — A/") In ( 1 — Af) , 
with J\f < 1 the mean occupation of the cavity. As usual, the 0D cavity is understood as a cavity 
of arbitrary geometry which can accommodate one particle at the most. 

To determine Kij(r) we will follow another procedure: we will impose that the low density 
limit of the second functional derivative of © with respect to the density profiles Pi(ri) and Pj(v 2 ) 
coincides with the overlap function of two HD of radii Ri and Rj, which turns out to be the exact 
low density limit of minus the direct correlation function, i.e. 

- r 12 ) = I dr' (^ (0) (r>f V12 - r')) + 2Jf y (r 12 ) / dv'uf \v')uf ) (r 12 - r'), (8) 

where (fy) = f i3 - + fa has been introduced to denote symmetrization of fa with respect to its 
indices, and R^' = Ri + Rj. The calculation of the integrals involved in Eq. ([8]) leads to 



K l3 (r 12 ) = Trr^sin- 1 t l3 \r l2 ))R^ I - t^f Q (r 12 - R^ 6 - r 12 ) , (9) 

, x A2 + Rl - R 2 i 

t,(r 12 ) = 12 2ri ^ >- , (10) 



where r\- ' = \R — Rj\. The kernel © is symmetric with respect to the exchange of indices 



and j due to the equality hi* = Ri^l — ^(r^) 2 = Rty 1 — £#(7*12) 2 » which is easily visualized 



FIG. 1: Sketch of the triangular geometry defined by the lengths Ri, Rj and r\ 2 . 



in Fig. Q] this figure shows a sketch of a typical configuration of two HD with different radii for 
which the Kij{ri 2 ) is different from zero; in it the height hij of the triangle formed by the lengths 
ri2, Ri and Rj can be calculated either as Ri sin0j or as Rj sincfij, thus proving the symmetry 
Kij(ri 2 ) = Kji(ri 2 )- Also, from Eq. © and the triangular geometry the kernel can be rewritten as 
Kij(ri 2 ) = 7TAij(r 12 )(j)ij(ri 2 ), with Ay = r^h^ the sum of the areas of both triangles (the shaded 
region of Fig. [B, and fcj = (sin -1 1^- ) , the angle formed by the sides Ri and Rj. 
For one component, Kij(r) recovers the expression 



reported in Ref. It is interesting to note that, because of the presence of weights u\ (r) 

(proportional to Dirac delta functions) in the definition of N(r), its expression can be greatly 
simplified. After insertion of (r) into Eq. © and integration over the radial variables n and r 2 
we find 

N(r) = — V RiRj / d<Pi / # 2 pi(r + RiU^p^r + i? j u 2 )T(0 12 ), (12) 
47r i,j J o Jo 

where 12 = (pi — 4> 2 , Uj = (cos sin (pi), and the function T((p) is defined as 

T(0) = |0-27rn||sin0|, (13) 

with n the integral part of the fraction (0 + 7r)/(27r). The first factor on the right hand side of 
(fT3l is the 2tx -periodic function shown in Fig. |2] with a dashed line. The function T(cp) is plotted 



.5 1 1.5 2 

<j)/(27t) 



FIG. 2: The function T(<p) in the interval [0, 4-7r] (solid line). Also is shown with dashed line the triangular 
wave function. 

in the same figure. From the new form of N(r) given by Eq. (TT2]) we can conclude that, despite 
the presence of a two-particle weight in its definition, the numerical cost required to evaluate it is 
the same as that required to calculate the local packing fraction 712(1"), because both quantities are 
defined through a double integral. 

Taking into account the result f dcb-\ f dcboT (Sio) = An 2 we find that the uniform limit of $( 2D ) 

: J J nnnn 

[cf. Eq. (OQ)] coincides with the scaled particle theory (SPT) result for a mixture of HD Illl.l2l.l3l.ll4ll 

$( 2D ) = -eoln(l-6) + ^T^V 5 (14) 

47T 1 — £2 

where £ = E* Pi, 6 = E* PifaRi), and 6 = E* Pi (^f )• 

Because the derivation of -fQj (r) has not followed the requirement of exact reduction to OD 
cavities, as in Ref. Jll| . the question arises as to whether this nice property holds for this new func- 
tional. In Appendix IA1 we proof a stronger property: the functional (OQ) fulfills an exact 2D— >1D 
crossover; in other words, by inserting the profile Pi(r) = pi(x)S(y) into ® and © we recover the 
exact functional for ID hard rod mixtures |3^], cf. Eqs. (IA16I) and (IA18I) . With this we have also 
proven that the present method and the cavity formalism of Ref. [11] are two equivalent methods 
to obtain a FMF for mixtures of HD. 
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III. DENSITY FUNCTIONAL FOR A MIXTURE OF PARALLEL HARD CYLINDERS 



In this Section we will construct a FMF for a mixture of PHCL starting from two different 
density functionals for HD mixtures: the first one is calculated through Eq. (0Q), while the second 
one is Rosenfeld's proposal yO 



$g D) = -n ln(l - n 2 ) + —- , (15) 

47r 1 — n 2 

where v = J2iPi* w i( r ) an d v = J2i Pi * w i( r )' with me new weights 



Wi{r) = 27cRiUjf\r), W;(r) =w~. (16) 

Note that the weighted density v (r) is not the absolute value of v(r). Rosenfeld obtained this 
expression by approximating the Mayer function of two HD by the sum of convolutions between 
single-particle weights fij(r) = (oof^ * u>y)(r) + (2ir)~ 1 [w i * Wj(r) — w, ; * Wj(r)] and requiring 
also that the scaled particle result (TT4l) was recovered in the uniform limit [3]. The expression (fT5l) . 
however, does not conform with any dimensional crossover to lower dimensions. 

The FMF for PHCL is obtained by resorting to the dimensional crossover 3D^2D which any 



functional (3f£ D) [{pi}} fulfill. In Ref. tl8|] it was argued that having a FMF for a mixture of two- 
dimensional particles JF e x D ^ [{pi}] one can construct an explicit expression for a FMF for a mixture 
of parallel anisotropic three-dimensional (3D) bodies whose constant section perpendicular to their 
main axes is that of the 2D particles. In the same reference it is shown that the resulting functional 
fulfills by construction the 3D^2D dimensional crossover when the original fluid is confined such 
that the centers of mass of the particles are confined on a plane perpendicular to their axes. We 
will apply this method to obtain a FMF for a mixture of parallel cylinders as follows. First of 
all we need to redefine the functional in such a way as to include the 2-coordinate dependence of 
the density profiles and correspondingly of the weighted densities. The new weights are obtained 
multiplying the old ones by the factors 0(Lj/2 — |^|), i.e. 

uf \v) = ^ (t V)0(V2 - |*|), ^ (3) (r) = V)0(V2 - |*|), (17) 
^\t u t 2 ) =n ij {ri, ri)Q(Li/2- 1^1)6(^/2- \z 2 \). (18) 

The resulting free energy density, which we denote $( 2D ), is the same as that given by Eq. (Q~|), 
but with the substitutions n (r) — > ni(r), n 2 (r) — > n 3 (r) and N(r) — > N 2 (r). The new weighted 
densities are obtained through the same expressions given for the 2D case but using the new 
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weights (fT71) . (fT8l) . Note that the vector position r is now defined as r = (i- 1 , z). The dimensional 
crossover 3D— >2D holds if the 3D excess free-energy density is obtained by the following formula 



y 



* (3D, (r) = £^:* <2D, M' (19) 

i 

In our case this leads to 

$< 3D ) = -n In 1 - n 3 + + 2 (20) 

l-n 3 (l-^ 3 ) 2 

where the one-particle weighted densities ri a (r) are calculated as usual as n a (r) = J2i Pi * ( r ) 
[* stands for the convolution f *g(r) = J dr' f(r')g(r — r')]. The one-particle weights are defined 
jointly by Eq. (TT7T) and 

4V) = ^ (0) (ri)i(V2 - |*|), 4V) = \4 \r ± )6(L l /2 - \z\). (21) 
The two-particle weighted densities are calculated as 

N «( T )=JL j rfr i / dv 2Pi {r 1 ) Pj {r 2 )^ ) {v-v u v-v 2 ), (22) 

with Q\f (r 1( r 2 ) given by Eq. (TT8l) while U-^ (ri, r 2 ) is 

ngW 2 ) = nvtfrf) - \ Zl \)5{L 3 /2 - |* 2 |)\ . (23) 

As mentioned above, the free-energy functional JF( 3D ) has a correct dimensional reduction 
to JF( 2D ), the free energy for a HD mixture, when the density profiles are taken as Pi(r) = 
p( 2D ) ^ r _i_^ _ (x ; y)], i.e. projecting the cylinders on the plane perpendicular to their 

axes. It is easy to show that the dimensional cross-over 3D— >2D, where the projection is now 
in a plane parallel to the cylinder axes, also holds. To show this we take the density profiles as 
Pi(r) = pi(r^)8(x) [r'l = (y, z)], with Pi(r-H) = pf D \y,z). Inserting these Pi (r) into $( 2D )(r) 
and using the already shown dimensional cross-over 2D— >1D of a FMF for a mixture of HD (see 
Appendix |Al). we obtain, from Eq. (TT9b . 



the free-energy 
their axes) 118 



$ (3D) (r) - *ffl(*) = -no(r) In [1 - n 2 (r)] + "f^^ , (24) 

density of a mixture of parallel hard rectangles (the section of the cylinders along 
, 19(1 . The weighted densities for such particles are now defined as n (r) = 



FIG. 3: Sketch of all dimensional crossovers fulfilled by the fundamental-measure functional for a mixture 
of parallel hard cylinders: (a) from cylinders (3D) to rectangles (2D), (b) from cylinders (3D) to disks (2D) 
and (c) from disks (2D) to rods (ID). 



E,pf D) * MO = E,pf D) * n lx {v) = £ 4 pf D) * [k t (y)M*)l. 

andni 2/ (r) = £ i pf D ^ * [9 ai (y)9L i {z)}, where the shorthand notations = |5(wi/2 — \s\) and 
^ui(s) = @(ui/2 — \s\), with s = y, z and ex, = 2i? i5 have been used. 

Thus the FMF for PHCL that we have just obtained conforms with all dimensional crossovers 
to lower dimensions, which we have sketched if Fig.O Nevertheless, this functional is not perfect 
because it shares with that of HD the defect caused by the existence of "lost cases" [U, lfl. There 
are three-point OD cavities such that particles sited at those points have pairwise overlap but no 
triple overlap. For those cavities the FMF of a HD mixture does not reduce adequately to $(° D ) (the 
lost cases for the one-component fluid were already pointed out in Ref. [1 1]). As a consequence 
of that the FMF for a mixture of PHCL suffers from the same illness. 

The excess free-energy density of PHCL obtained from its HD counterpart using (TT9b and 
Rosenfeld' approximation $^ D ' ) [cf. Eq. (fT5l) l results in 



.(3D) W1 v . n x n 2 + v lV 2 - viv 2 1 n 2 (vj - v 2 2/ 

= ~ n o In 1 - n 3 ) + + — , (25) 

l-Tia 4vr (1 - n 3 ) 2 

where v a (r) = Y^iPi * wf\r), v«(r) = J2 i p l * wf } (r), and 



Wi (1) (r) = w 2 (1) (r) = ^\rf-, (26) 



R 

wf\v) = 2nR4 1 \r), wf } (r) = ^(vf-, (27) 

rii 



where the uo^f 1 (r) are those defined by Eqs. (flTT) and (T2T|) . 
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IV. UNIFORM MIXTURES 



In this section we give the explicit expression for the uniform limit of the FMF for a mixture of 
PHCL, which coincides with the SPT result. 

It is easy to show that J dri J cfa^f^., (ri, r 2 ) = nRiRj. Taking into account this result, the uni- 
form limit jOj(r) = pi of both free-energy densities, $( 3D ) from Eq. (l20l) and $^°' ) from Eq. (|25|) . 
yield the result 



< 3D) = -eo in(i - 6) + + rjr^ 

1-43 87T (1 - 4 3 ) 



(28) 



where we have defined the vectors £ t = £|j (i — 1, 2), with components 

^ = ^>^, d = X>y' (29) 

^ = ^ piiixRiLi, d = J2 p* 2nR ^ (30) 

i i 

while £ = Z)j P« an ^ £3 — Si Pi^RjLi are the total density and total packing fraction of the 
mixture, respectively. From Eqs. (f28l) — (f30b we can see that the excess part of free-energy density 
is a function of certain weighted densities Q a \ which can be calculated as the sum of products 
between the particle densities pi and their fundamental measures: {Ri, Li/2}, the principal radii 
in the directions perpendicular and parallel to the cylinder axes, {2nRiLi, 2nR}}, the areas of the 
surfaces oriented along the perpendicular and parallel directions, and nR^Li the particle volume. 

Within the SPT formalism, the excess part of the free-energy density of any mixture of convex 
particles should fulfill the following differential equation Jj, LL90 



$ + 2^^— + £o = — • (3D 



This equation holds for (1281) . thus showing that of our functional gives the SPT result. Finally, the 
equation of state within SPT can be calculated as (3P = -^7—, resulting in 

l-£ 3 (l-£ 3 ) 2 47r(l-£ 3 ) 3 
This equation of state can be used to study the possible demixing scenarios that a mixture of PHCL 
has. 

In order to show the existence of demixing in a binary mixture of PHCL we first specialize the 
mixture to the case in which both particle volumes are unity, i.e. vx — v 2 — 1. This assumption 
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allows us to calculate the particle lengths and diameters ■ ' and <jj = ck- 1 , where 

Ki = Lijoi are the cylinders aspect ratios and c = (4/vr) 1 / 3 . It is easy to show that the expression 
(|28T) for this particular mixture gives us the following expression for the free-energy per particle 

<p = ($ id + $1 3D) ) jp (with $ id = ]r\ Pi IHviPi) - 1]) 

(p = lay — 1 + x\nx + (1 — x) ln(l — x) + yS(x; r) + y 2 T(x; r), (33) 
while the expression for the fluid pressure is 

(3Pv i = y + y 2 S(x;r) + 2y 3 T(x;r) } (34) 

where y = 77/(1 — 77), 2 = X2 is the molar fraction of species 2 and r = k 2 /ki is the ratio between 
the particles aspect ratios. Note that for this particular mixture we have 77 = p. Also, 



S(x; r) = 3 + (r 1/3 - l) (l + 4r" 1/3 + r" 2/3 ) z(l - x), (35) 
T(ar; r) = 1 + (r 1/3 - l) 2 [r" 1/3 (2 + r" 1/3 ) + (l - r~ 2/3 ) x] x(l - x). (36) 



Note that while the function S(x;r) = 5(1 — x; r) is symmetric with respect to the value x — 1/2, 
T(x; r) is not. Thus the spinodal instability curve with respect to the phase separation is not 
symmetric with respect to x = 1/2. Besides we have the obvious symmetry S(x;r) = S(l — 
x; r" 1 ) and T(x; r) = T(l — x; r _1 ). 

The lost of mixture stability with respect to phase segregation can be calculated as usual as 

( «9 2 $ (3D) \ 

del + w)= ' (37) 

which is equivalent to the following condition, expressed in the variables y and x, 



d_ 

dy 



y 2 ir 



d ^ fy^) =0. (38) 



dx 2 \ dydx / 

Inserting (l33l) into ((381) we calculate the demixing spinodals for different values of the asymmetry 
parameter r in the plane x — 77. Figure H] shows these demixing spinodals for values r = 20, 10, 5 
and 2. 

Of course this analysis does not prove that a thermodynamically stable fluid-fluid demixing 
occurs, as inhomogeneous phases are not being accounted for. 
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FIG. 4: Demixing spinodals for the phase separation between two nematic phases of different composition 
calculated for values of the asymmetry parameter r as labeled in the figure . 



V. DIRECT CORRELATION FUNCTION 

The second functional derivative of ftFe^ = J dr $( 3D )(r) with respect to the density profiles 
Pi(r) and Pj(r) evaluated at the uniform densities pi gives us, after a long and tedious calculation, 
the following expression for the direct correlation function 



where Xi 



Cii(ria) = [Xo + Xi • AR^r^) + % 2 • AS i3 -(r 12 ) + XaAV^ru)] fij(r 12 ) 
d{pP) 



(39) 



and 



h(ru) = -6 - r 12 ) 6 (4 +) /2 - \z 12 



AR ?; ,(r 12 ) 
ASy(r 12 ) 
A^-(r 12 ) 



ALj j (r^)/(2n),A4 3 (z 12 )/2 



A^(r 12 ),A4(^) 
-A4(ri)A4(^ 12 ) 



(40) 
(41) 
(42) 

(43) 



with L^f = \Li±L j\, are the Mayer function (1401) and the geometrical measures of the body 
defined by the overlap between two cylinders i and j whose centers of mass are separated by the 
vector r 12 . These measures are characteristic radii along the perpendicular and parallel directions 
(|4"T1) . the oriented surfaces (l42l) and the total overlap volume (|43T) . The radii in turn are defined 
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through the total length of the arches and the height of the overlap body. The expressions 
for these quantities as well as for AS^ (a =_L, ||) are given in Appendix [Bj 

The form of the direct correlation function (l39l) as a function of the geometric measures of the 
overlap body is exactly the same as that obtained from the Percus-Yevick approximation for a HS 
mixture, as it was first shown by Rosenfeld 3]. The same formal expression is also obtained 
for a mixture of parallel hard cubes fl8[[l9| . 



VI. CONCLUSIONS 



We have derived a FMF for a mixture of HD, and further used it to construct another one 
for a mixture of PHCL. The resulting functional fulfills all dimensional crossovers, a feature that 
makes the obtained functional very useful in the study of fluid mixtures of perfectly aligned hard 
rods confined by external potentials. These external potentials may have planar or cylindrical 
geometry. Some interfacial phase transitions, such as wetting, layering and capillary ordering, can 
be studied as well using this functional. 

Of course, the parallel alignment constraint limits the use of the PHCL functional to the study 
of highly oriented phases, such as nematic, smectic or crystal phases at very high pressures. A 
particularly interesting application of this functional is the determination of the phase behavior of 
polydisperse hard rod mixtures. The inclusion of smectic and columnar phases in the study makes 
the constraint of perfect particle alignment indispensable to achieve the numerical minimization 



of the functional. Some experimental works [|37L 13 81] as well as simulations Q40] predict that 
polydispersity enhances the columnar phase stability with respect to the smectic phase. It will be 
interesting to check these conclusions with the reported functional. 
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APPENDIX A: 2D^1D LIMIT OF THE FMF FOR A MIXTURE OF HARD DISKS 



We begin with the calculation of the one-dimensional limit for the two-particle weighted density 
N(r) defined in Eq. ©. Substituting the expressions p,(r) = pi(x)5(y) (where Pi(x) is the one- 
dimensional density of species i) and integrating over the coordinates yi (i — 1, 2) we obtain 

N ( x > y) = ^~2 ^Tr - / dxi / dX2 Pi( x ^Pj( x ^) 6 [ R i - V ( x - x i) 2 + y 2 ) 

x S (Rj - y/(x- x 2 ) 2 + y 2 ^j K i:j (\xi - x 2 \) 



y y j ^— ^ — y^^i — yj. i i dx 2 Pi{xi)pj(x 2 ) 

j ■ UiUj J-oo J-oo 

x 5{ui — \x — xi\)S(uj — \x — x 2 \)Kij{\xi — x 2 \), 



(Al) 



where Ui = a/ R 2 — y 2 and we have used the identity 

6 (Ri - ^(x- Xl ) 2 + y 2 ) = &(R t - \y\)—5( Ui - \x - Xl \). (A2) 
Because of the deltas in the integral (IA1I) . it can be readily performed and yields 

N(x, y) = S~] — H) e R 3 V_ _|_ Ui )pj(x + uj) + pt{x - Ui)pj(x - Uj)]Kij(\ui - i 

i,j J 

+ [pi(x + Ui)pj(x - Uj) + Pi{x - Ui)pj(x + Uj)]Kij{Ui + Uj)} . 

(A3) 

From Eqs. © and (PT0l) for K^r) we obtain 

i^d^i ± = 7r|y||ui ± Uj\\Xi ± Aj|, (A4) 

where we have defined 

Xi = sin-^Ui/Ri) = coe-Wvl/Ri). (A5) 

In order to proceed let us assume for a while that Ri > Rj. Then Ui > Uj and Aj > Aj and 
therefore 

9(^-1^1)0(^-1^1)1^/1 



N(x,y) =J2 



Air 



x | [pi(x + + pi(x - Ui)] [pj(x + Uj) + pj(x - uj)] (J^- + ^-J (A6) 

- [pi(x + Ui) - Pi(x - Ui)][pj(x + Uj) - Pj(x - Uj)] (— + ■' 



Ui Uj 
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In order to get the equivalent expression when R { < Rj we should just exchange the indices i and 
j in the above expression. But this expression is invariant under this exchange of indices, therefore 
it holds for any Ri and Rj. 

Let us now obtain the densities n (r) and n 2 (r) given by Eqs. © and ©. When inserting the 
one-dimensional density profile one gets 



n 2 (x,y) = ^n 2i (x,y), 

i 

n o(x,y) = ^2n 0i (x,y), 



n 2i (x,y) = Q(Ri - \y\) 



X+Ui 



Pi(t) dt, 



(A7) 



n 0i (x,y) = — [pi{x + Ui) + pi(x - (A8) 

ZlTUi 



Notice that, as = Ri when y — 0, then n 2 (x, 0) = n\{x) and n 0i (x, 0) = n 0i (x)/iiRi, with 



x+Ri 



piit) dt, 



n (x 



x—Ri 



noi(x) = - J2iP*( x + R i) + Pi( x - ( A9 ) 



the two weighted densities of the exact DF for a mixture of ID hard rods [39]. 

For the sake of notational clarity, in what follows we will omit the arguments of n (x, y) and 
n 2 [x, y). Equations (IA7I) and (|A8I) help us to rewrite (IA6I) as 



N(x,y) = Y \ 7r\y\n 0i n oj (uiXi + ujXj) 

i,3 



\y\ dn 2i dn 2j ( A; Xj 



Att dx dx V Uj u 



Ely I dn 2 A,' dn 
UiXinoi - — ^- 2^ 



(A10) 



27T <9lC ^— ' Mi <9x 



(for notational simplicity we have omitted the x and y dependence of the weighted densities). 
Now we can integrate $^ 2D )(x, y), as given by Eq. ©, with respect to y to obtain 

N(x,y) 



/OO /»C 
& m \x,y)dy = / 
-oo J — c 



dy < -n (x, y) ln[l - n 2 (x, y)] + 



1 - n 2 {x,y) 



iV> )2 /) 



(All) 



n 2 {x,y) 

because the integrand is an even function of y. On the other hand, when y > 

<9n 2 



9j/ 



-27ryn , 



5 A^y <9 2 n 2i 
— (MiAjnoj) = -not - - ^-5- 



(A12) 
(A13) 



therefore 

$( 1D )(x) 



dy 



dyy^AlT. K-^oi ln(l - n 2 ) 



<9x 



-^-lnfl - n ) dU2i 
2nUj dx 



(A14) 
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The first term in this equation can readily integrated. Since for y = we have «i(0) = Ri and 

Aj(0) = 7r/2 and for y = Ri we have Ui(Ri) = Aj(-Rj) = 0, it follows that 



where 



ox 



$( 1D )(x) = -n {x) ln[l - m(x)} 



(A15) 



(A16) 



is the exact DF for a ID hard rod mixture [[390 in terms of the weighted densities (|A9I) . and 

dn 2i 



poo 
./0 



m(l - n 2 )—— 

TTUi OX 



(A17) 



Assuming proper boundary conditions for the density when x — > ±oo, the free-energy functional 
for the system is given by 



/3^L 2D) [{A}] 



¥ 1D \x)dx= / & 1D \x). 



(A18) 



This completes the proof of the exact 2D^1D dimensional crossover of the DF for HD (OQ). 



APPENDIX B: GEOMETRIC MEASURES OF THE OVERLAP BETWEEN TWO CYLIN- 
DERS 

In this appendix we provide explicit expressions for the geometrical measures of the body 
formed by overlapping two cylinders with radii Ri and lengths L^. To begin with, the formal 
definition of all these measures is 

&(ria) = -6 R ( + )(ri 2 )Q L(+) (z 12 ), (Bl) 

A4(r 12 ) = Q RW (ri 2 ) [Q Li * e Lj (z 12 )] , (B2) 

ALJ-Cna) = (6* * 5 R .(ri 2 )) 6 l(+) (z 12 ), (B3) 

A4(r 12 ) = 2 [6* * 0a, (r£)] (B4) 

where 6 u (r ± ) = Q(u — r- 1 ), u (z) = Q(u/2 — z), and 5 u {r L ) = 5{u — r L ). It is rather easy to 
evaluate these expressions appealing to their geometrical meaning. Hence, the total arch lengths 
of the cross-section of the overlap body is 

AL^) = 2 { (Ri cos' 1 ty) (r£ - i^) + | (^ +) - i£>) 6 (4") - r^ 2 ) } , (B5) 
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while its height is 



a4(* 12 ) = 



-(+) 



L 



(-) 



\ Z 12\ 



z 12 \\q{l^/2-\z 12 \). 



2 1 I 2 

Similarly, the expression for twice the area of the base of the overlap body is given by 



a4(^ 2 ) = 2 (if: 



cos 1 tjj — 



12 

2 



+iK'-^)' e (^-' 1 



12 r> 



while its lateral area is readily obtained as 



(B6) 



(B7) 



A^(r 12 ) = ALi(^ 2 )ALf.(z 12 ). 



(B8) 



[1] Y. Rosenfeld, Phys. Rev. A 33, 2025 (1986). 
[2] Y. Rosenfeld, J. Chem. Phys. 89, 4272 (1988). 
[3] Y. Rosenfeld, Phys. Rev. Lett. 63, 980 (1989). 
[4] Y. Rosenfeld, Phys. Rev. A 42, 5978 (1990). 
[5] Y. Rosenfeld, J. Chem. Phys. 98, 8126 (1992). 

[6] A. Gonzalez, J. A. White, F. L. Roman, and R. Evans, J. Chem. Phys. 109, 3637 (1998). 
[7] A. Gonzalez, J. A. White, F. L. Roman, and S. Velasco, J. Chem. Phys. 125, 064703 (2006). 
[8] P. Tarazona, Phys. Rev. Lett. 84, 694 (2000). 

[9] Y. Rosenfeld, M. Schmidt, H. Lowen, and P. Tarazona, J. Phys.: Condens. Matter 8, L577 (1996). 
[10] Y. Rosenfeld, M. Schmidt, H. Lowen, and P. Tarazona, Phys. Rev. E 55, 4245 (1997). 
[11] P. Tarazona and Y. Rosenfeld, Phys. Rev. E 55, R4873 (1997). 
[12] P. Tarazona, Physica A 306, 243 (2002). 

[13] R. Roth, R. Evans, A. Lang, and G. Kahl, J. Phys.: Condens. Matter 14, 12063 (2002). 
[14] H. Reis, H. Frisch, and J. L. Lebowitz, J. Chem. Phys. 31, 369 (1959). 

[15] J. A. Cuesta, Y. Martinez-Raton, and P. Tarazona, J. Phys.: Condens. Matter 14, 11965 (2002). 
[16] J. A. Capitan and J. A. Cuesta, Phys. Rev. E 76, 011403 (2007). 
[17] J. A. Cuesta, Phys. Rev. Lett. 76, 3742 (1996). 

[18] J. A. Cuesta and Y. Martinez-Raton, Phys. Rev. Lett. 78, 3681 (1997). 

18 



[19] J. A. Cuesta and Y. Martinez-Raton, J. Chem. Phys. 107, 6379 (1997). 

[20] L. Lafuente and J. A. Cuesta, J. Chem. Phys. 119, 10832 (2003). 

[21] L. Lafuente and J. A. Cuesta, Phys. Rev. E 68, 066120 (2003). 

[22] L. Lafuente and J. A. Cuesta, J. Phys.: Condens. Matter 14, 12079 (2002). 

[23] L. Lafuente and J. A. Cuesta, Phys. Rev. Lett. 93, 130603 (2004). 

[24] L. Lafuente and J. A. Cuesta, J. Phys. A: Math. Gen. 38, 7461 (2005). 

[25] Y. Rosenfeld, Phys. Rev. E 50, R3318 (1994). 

[26] A. Chamoux and A. Perera, J. Chem. Phys. 104, 1493 (1996). 

[27] M. Schmidt, Phys. Rev. E 63, 050201(R) (2001). 

[28] J. M. Brader, A. Esztermann, and M. Schmidt, Phys. Rev. E 66, 031401 (2002). 

[29] A. Esztermann and M. Schmidt, Phys. Rev. E 70, 022501 (2004). 

[30] A. Esztermann, H. Reich, and M. Schmidt, Phys. Rev. E 73, 011409 (2006). 

[31] Y. Martinez-Raton, Phys. Rev. E 69, 061712 (2004). 

[32] Y. Martinez-Raton and J. A. Cuesta, Phys. Rev. Lett. 89, 185701 (2002). 

[33] Y. Martinez-Raton, J. Chem. Phys. 118, 10164 (2003). 

[34] M. Bier, L. Harnau, and S. Dietrich, Phys. Rev. E 69, 021506 (2004). 

[35] Y. Martinez-Raton, Phys. Rev. E 75, 051708 (2007). 

[36] L. Harnau and S. Dietrich, in Soft Matter, edited by G. Gompper and M. Schick (Wiley- VCH, Berlin, 
2007), p. 159. 

[37] F. M. van der Kooij, K. Kassapidou, and H. N. W. Lekkerkerker, Nature 406, 868 (2000). 
[38] F. M. van der Kooij and H. N. W. Lekkerkerker, Phys. Rev. Lett. 84, 781 (2000). 
[39] T. K. Vanderlick, H. T. Davis, and J. K. Percus, J. Chem. Phys. 91, 7136 (1989). 
[40] J. M. Poison and D. Frenkel, Phys. Rev. E 56, R6260 (1997). 



19 



